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ABSTRACT 

We have constructed a Hadamard Transform Spectrometer 
(HTS) which simultaneously obtains fifteen infrared spectra, 
each having 255 spectral elements. Spectra are obtained 
essentially in real time through use of a minicomputer with 
8 k words of memory and a CRT display. This permits opera- 
tion of the instrument in the field. 
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1. Introduction 

The principle and operation of Hadamard Transform 

Spectrometers (HTS) has by now been described in a number 
12 3 4 

of papers * . The theoretical predicted multiplex ad- 

vantage has been experimentally verified for one of the 
configurations in which these instruments can be operated , 

r -r 

and astronomical observations have been carried out° . 

There are basically three modes in which these instru- 

12 3 

ments can operate. The simplest model * ’ has a single 

entrance slit and n exit slits. It produces a spectrum 

3 4 

having n spectral elements. A more complex instrument-'* 
is one having m entrance slits and n exit slits, m < n. 

It can be operated in two different modes. If the radiation 
incident on the spectrometer is diffuse and homogeneous, the 
prime advantage of having several entrance slits is an in- 
crease in the instrument's throughput j but the final pro- 
duct still is a single spectrum having ~ n spectral elements. 
Alternately, if the radiation incident on the instrument 
is not homogeneous, one may wish to obtain a separate spec- 
trum for light incident on each entrance slit position. 

One application for such an instrument that we hope to 
realize is a simultaneous analysis of the spectra of the 
different cloud bands of Jupiter. By aligning the instru- 
ment's entrance slits parallel to the Jovian equator, any 
marked differences in the infrared spectra at different 
Jovian latitudes should become apparent. 

This last mode of operation is the one viiich we have 
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been able to make practical, in the sense that m and n 
are quite large, that the equipment is moveable, and that 
the operation of the instrument is essentially instantaneous. 

The only previously discussed instrument of this kind 

O 

was described by Phillips and Harwit , who constructed a 
pilot model having .19 exit slits and up to 19 entrance 
slits. With this instrument these authors were able to 
demonstrate the principle of operation. 

The new instrument we describe has 15 entrance slits 
and 255 exit slits. It obtains 15 different spectra each 
having 255- spectral elements. The total number of intensity 
measurements that have to be made is 15 x 255 = 3825. Data 
can be gathered at a rate of up to five readings per second 
so that the entire spectral run lasts only some thirteen 
minutes. This rate - , which is perfectly adequate for most 
astronomical observations, is limited primarily by the rate 
at which our minicomputer can process the incoming data 
in real time. At the end of the thirteen minute run, the 
spectrum is available for immediate viewing on a CRT display. 
In our laboratory we photograph this display with Polaroid 
film, but hard copies of other types could also be used. 

In addition we can punch and store the spectra on paper 
tape for further processing. 

2. Instrumentation 

The optical configuration, the drive of the encoding 
masks that are used and the spectral capabilities of our 
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instrument have been described in detail by Harwit, et al. , 
in sections VII A and B of their paper. For brevity we will 
not repeat this description, but we are including a schematic 
diagram of the optical layout. Fig. 1., and a table of 
instrumental parameters. Table I. 

The actual operation of the instrument procedes in 
exactly the manner specified there. The only difference 
lies in the way that the data points are handled. Harwit, 
et al. , used these data to obtain a single spectrum averaged 
over the entire radiation collected by the fifteen entrance 
slits of the instrument. Here we describe the handling of 
the data to obtain fifteen different spectra, one corresponding 
to each entrance slit of the spectrometer. 

The main difficulties that arise in the implementation 
of this mode of operation involve speed, ease of verification, 
and self-containment. We wanted to construct an instrument 
that could be transported to any astronomical observatory 
where infrared observations can be made, and we wanted to 
be able to verify in real time that adequate spectra were 

f 

being obtained. We therefore chose to center our construc- 
tion around a readily transported minicomputer programmed to 
produce a cathode ray display. 

We have used this display primarily, for rapid visual 
inspection and for photographs that can be used as a quick 
reminder about information punched out and stored on paper 
tape. Since punching 5825 spectral points on paper tape 
requires some 1 J> minutes, we plan to punch and store only 
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those spectra that are of genuinely high quality to warrant 
further processing in the laboratory. Processing of this 
type may involve various atmospheric corrections, normali- 
zation to a calibration standard, and so on, 

J. Algorithm for Real Time Decoding of Spectral Data 

In our minicomputer assisted system we make use of the 

following algorithm for reducing the spectra in real time. 

Let both the entrance and the exit masks be linear 

arrays encoded by Reed Muller codes'*'. Then the matrix 

of spatial-spectral elements f is related to the matrix of 

-1 

measurements t\, by and postmultiplying by S 

s^S = tj (1) 

To obtain the spatial-spectral information about the 
viewed scene we solve this equation by premultiplying the 
data matrix by s~"*" 

ij/ — s "*" . (2) 

Now consider the element r^. It is multiplied only 
by elements of the first column of s~" L ; and in turn it 
multiplies only the elements of the first row of S \ To 

a given spectral-spatial element it therefore contributes 

-1 -1 -1 -1 
an amount s^ • But the elements s^ and S-^j 

all have values, either +1 or -1, and the result is that 

each element ■$. . of the matrix i/ receives a contribution 

+T)n or -t) u from the reading r^. 

This procedure is generally valid. Any reading 

will make additive -contributions that can only have values 
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+ r) k ^ or -r) k ^, to each element of the y matrix. 

For real time decoding we therefore need the fol- 
lowing: 

(a) A memory that consists of tins containing the 
contributions to the elements ijf^^ accumulated up to any 
given time, t, in the cycle of measurements. For a device 
that can resolve m spatial and n spectral elements, this 
memory requires mn bins and of the order of mn memory words, 

(b) For each acquired reading we perform a series 
of additions of values either +T) k ^ or -r^, one to each of 
the mn memory bins f. .. But before that can be done, we 
need to decide on the assignment of + or - needed for a 
given bin. This is done in the following way: 

We store the sequence of + and - signs in one column, 
of s ^ and in one row of S ^ and in one cycled permutation 


of each of these vectors. Let us designate these signs 

by their positions in these two vectors, as s ." 1 and S . -1 

J 

respectively, i=l,...mj J=l,.,.n. (Since each of these 


sequences is cyclic it can respectively be brought into 
its and cycling position after a measurement r^. ) 
The elements of the two vectors then are multiplied in all 
possible combinations, to give a matrix having mn elements. 


y.. = s.-v - 1 
• 1 0 


i ’ 1 f # » # tn. i j ~ * 1 jf » * • n 4 


Each element . is either + or - depending only on 

whether the signs s . -1 and S . -1 are similar or dissimilar 

^ J 

for a particular combination of i and j values. 
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The additions h-t]^ or -tj^ to the bins are made 

as successive elements ) . . are computed, so that the elements 

L-i±J 

j need never be stored. Fig. 3 shows the relation of 
Y a j to a super-array containing the set of all elements 
that are constructed at various stages of. the computation.. 

When only a restricted number of spectral elements 
is of interest, we need to compute elements ^ . . representing 
only selected j values. This might be useful, for example, 
if only certain atmospheric COg lines needed to be -studied, 
and the spectral elements inbetween were of lesser interest. 
In that case only elements corresponding to a given 

j values need to be used, and the computing time decreases 
as p/n, where n is the total number of available spectral 
elements, and p is the actual number' of interest. 


4. Actual Use of the Computer 

We have used Computer Automation minicomputer models 
L. S.I. and also Alpha-l6 with the program we are about to 
describe. 

We start out with the inverse of the codes, s - ^ and 
S , for the entrance and exit masks stored in the computer. 
We store 509 elements of the exit mask, i.e. the 255 elements 
exhibited in the first mask position, plus an additional 
254 elements representing the further cycling of this mask. 
Similarly we store 29 elements for the entrance mask* 
representing the first 15 elements used, plus the further 
cycling of 14 elements. 
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For each reading we essentially make use of the 

matrix Fig. J>, making use of elements k to k+14 of the . 

* 

stored entrance code and elements £ to i+254 of the stored 

exit code. This matrix consists of + and - signs. When 

s ik _1 ^ ave the same sign, both being + or both 

being -, then the matrix position ij is assigned a + sign 

and the reading, is added to the accumulatively stored 

value of if/.. . If the elements s. and S,. -1 have dis- 

ij 

similar signs, then a - sign is assigned to ij and the 
reading is subtracted from the stored values. This 
whole process takes ~ 100 millisec, and is carried out while 
the succeeding intensity measurement r) k is being made. 

In practice we start with the first spatial element, 
i = 1, and add or subtract the contributions to all the ^ 

X J 

values, successively going from' j = 1 to j = 255. We then 
repeat this procedure for i values going from 2 to 15. 

This whole procedure is carried out while the exit mask is 
moving from position £ to £+1 depending on whether the 
exit mask is moving forward or back. The entire process 
is then repeated for the next reading . When the 

exit mask reaches its 255th position, £ remains unchanged, 
but the entrance mask moves from the position k to k+1. 

For odd values of k, the exit mask moves in the direc- 
tion of increasing £ values, and for even values of k, it 
moves toward decreasing values. In short, the exit mask 
moves back and forth as readings are taken. After the 
entrance mask has moved through all its 15 positions, and 
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the total 5825 readings have "been taken and added onto or 
subtracted from the elements, the run is completed. 

A display program allows us to display ( the information 
in a number of different ways: 

a) We can call for the spectrum corresponding to any 

one of the entrance slit position, and display it individually. 

b) We can display all fifteen of these spectra simul- 
taneously, with the zero baseline of each spectrum vertically 
displaced from the next one. While this format is somewhat 
crowded, it does allow a quick comparison of the individual 
spectra. 

c) Finally we can display the sum of the different 
spectra. In order to do this, we have to take into account 
that the spectrum for given entrance slit is displaced 

by one spectral element from the spectrum for the next 

O 

entrance slit position 0 . In other words, the wavelength 
for element corresponds to the wavelength for element 

j+l* because of the slightly displaced light paths 
through the spectrometer. 

These three forms of display are' shown in Figure 4 
as a), c) and b) respectively. 

We feel we have been successful in realizing an optical 
system that had been under discussion for some time^ 5 ^, 
and that the practical problems faced in constructing such 
a system are not excessive. 
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Table I. Relevant Instrumental Parameters 


Focal length 

= 49.5 cm 

Grating ruling = 20 lines/mm 

f number 

= 7.5 

Grating order = First 

Slot width 

= 0.1 mm 

Grating blaze - 9.5(d in first 

Slot height 

= 3.5 mm (or 

order 

Readings 
per second 

10 mm) 
= 5 

Range* = 8-l4|i (1250- 

•715 cm" 1 ) 

Resolution = O.Olp, (1 cm" 1 ) 


*The instrument can also be tested in the 1.7 |x region in 
6th order. 
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Figure Captions 

Fig. 1., Schematic layout of the Optical System 

Fig. 2. Explanation of the significance of the different - 
matrices and their rows and columns. 

Fig. 5» Matrices generated by the computer during the 
reduction of the spectral data. (See text for a 
full explanation. ) 

Fig. 4. Calibration spectra obtained for the mercury vapor 
emission lines at 1.6$ and 1.71p. (a) The eighth of 

a series of fifteen individual spectra obtained; 

(b) an average of all 15 spectra; (c) a display of the 
15 spectra, each spectrum being displaced vertically 
from the next. The diagonal pattern near. the right 
hand edge represents a displacement of the peak be- 
tweeen successive spectra. This represents the actual 
shift in spectral range between adjacent spatial 
elements . 




Post Optics 



Each column 
represents the 
code of one of the 
successive spectral 
encoding mask 



Each column represents 
the reading obtained in 
one of the successive 
spectral mask positions 



readings obtained in one 
of the successive spatial 
mask positions 


Each row represents a 
different spectral element 





2m — I 


column i 


Super- array of + and — signs. 

Row m + l to 2m— I are a repetition 
of rows I to m-l. Columns n + l 
to 2n-l are a repetition of 
columns I to n — I 


2rT^I 












